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I. INTRODUCTION 



The study of unbound or loosely bound exotic nuclei casts a new light on the discovery of 
nuclear halo structure [l| . The valence nucleons can be readily scattered into single-particle 
resonant states in the continuum. Due to the small binding energy and small or no centrifugal 
barrier conditions, the valence nucleons tunnel out of the potential barrier to long distances 
with an extended density tail to form the so-called halos. Therefore, the resonant states in 
the continuum and the coupling between bound states and the continuum near the threshold 

n rt 

play an important role in the description of halo phenomena |2|, |3[. In most calculations, 
the continuum is replaced by a set of positive energy states without the contribution of the 
widths, determined by solving the Hartree-Fock-Bogoliubov (HFB) or relativistic-Hartree- 

nn 

Bogoliubov (RHB) equations in coordinate space and with box boundary conditions |5[. 
The properties of the quasipaticle resonant states in the continuum, such as the widths, are 
taken into account for the neutron case by the resonant Bardeen-Cooper-Schrieffer approach 
based on the relativistic mean field theory (RMF-rBCS) js-S]. 

Neutron halos have been prominently observed in experiment loiji^l and investigated 

17| , while those for proton halos are still few because 



HI 



20|, especially for the two-proton hal o |21 |. As one of the 



for neutron rich nuclei theoretically 
of the Coulomb barrier 

candidates, .^£Ne aroused great interest from both the experimental [22 
aspects 



- |24l | and theoretical 

25|-|28|. Latest measurements of high-precision mass and charge radius on ^''"^^Ne, 



including the proton-halo candidate ^''Ne, have been performed with Penning trap mass 
spectrometry and collinear laser spectroscopy [23^. The charge radii of ^''Ne is determined 
for the first time, as well as the binding energy and root mean square (rms) experimental 
error bar. Recently, the reaction cross sections for ^''Ne are measured by transmission 
method and the corresponding density distribution is deduced from a modified Glauber- 



type calculation 



24|. 



In the present paper, we focus on the properties of the single proton resonant states, which 
are extracted from the analytical continuation of the coupling constant (ACCC) approach 
within the framework of the RMF theory. A theoretical description of two-proton halo can- 
didate ^''Ne is reasonably given, including the binding energy, matter radius, charge radius, 
density and occupation probabilities of the proton resonant states. The full self-consistent 
RMF-I-ACCC+BCS approach with the relativistic Coulomb wave functions boundary con- 
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dition are used and systematically introduced in Sec. II. The numerical details are given in 
section III. We apply this scheme to calculate the proton resonant states in ^^Ne and com- 
pare the theoretical results and the experimental results. Finally, we give a brief summary 
in Sec. IV. 



II. THEORETICAL FRAMEWORK 



A. The Relativistic Mean Field Theory 

The basic ansatz of the RMF theory is a Lagrangian density whereby nucleons are de- 
scribed as Dirac particles which interact via the exc hang e of various mesons (the scalar a, 
vector u and iso-vector vector p) and the photon 



29 



33| 



(1) 



where M is the nucleon mass and mg. (ga), fTi^ (gui), and nip (gp) are the masses 
(coupling constants) of the respective mesons. A nonlinear scalar self-interaction 
f/(cr) = -m^o"^ + ^cr^ + ^'^'^ meson has been included jsj]. The field 

tensors for the vector mesons are given as 

n^"' = d^'p" - d^p^" - gPip^" X p"), (2) 
The classical variation principle gives the following equations of motion 



[<x-p + Vv{r)+/3{M + Vs{r))]i:, = e,V^, 



(3) 



for the nucleon spinors, where Ei and ipi are the single-particle energy and spinor wave 
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function respectively, and the Klein-Gordon equations 

{-Aa + U'{a)) = g^ps, 

(-A + ml)uj^ = g^3^{v), 

(-A + m2) = gpf{r), 
-A <(r) 



(4) 



for the mesons, where 



Vv{r) = g^u'' + gpTspl + -eil-T3)A^, 



(5) 



Vs(r) 



are the vector and scalar potentials respectively and the source terms for the mesons are 



Ps 


A 
i=l 




A 
i=l 




A 

i=l 




i=l ^ 



(6) 



It should be noted that the contribution of negative energy states are neglected, i.e., the 
vacuum is not polarized. Moreover, the mean field approximation is carried out via replacing 
meson field operators in Eq. ([3]) by their expectation values, since the coupled equations 
Eq. (jS]) and Eq. (jl]) are nonlinear quantum field equations and their exact solutions are very 
complicated. In this way, the nucleons are assumed to move independently in the classical 
meson fields. The coupled equations can be solved self-consistently by iteration. 

For spherical nuclei, the potential of the nucleon and the sources of meson fields depend 
only on the radial coordinate r. The spinor is characterized by the angular momentum 
quantum numbers /, j, m, the isospin t = ±- for the neutron and the proton respectively, 
and the other quantum number i. The Dirac spinor has the form 



(7) 
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where Y-j^{6, 0) are the spinor spherical harmonics. The radial equation of the spinor, i.e. 



Eq. ([3]), can be reduced to 



(8) 



-|: + ^) Fl^ir) + [M + V,ir)] ^) = e^G^ir), 

(|: + ^) G^r) - [M - VUr)] F^^ {r) = b^fI^t), 

in which Vp{r) = Vv{r) + Vs{r), Kn(r) = W(r) - Vs{r), and k = {-iy+^+^/\j + 1/2). The 
meson field equations can be reduced to 

where (j) = a, u, p, and photon ( = for photon). The source terms read 

—QaPs — 920''^ (r) — g3(7^{r) for the a field, 



9pp3{r) 
epc{r) 



for the u field, 
for the p field, 
for the Coulomb field. 



and 



(10) 



(11) 



4vrrV.(r) = Ed^.W^ - |F,(r)p), 
1=1 

47rrV.(r) = E(|G,(r)p + |F,(r)n, 

i=l 

4vrrV3(r) = E (IG^WP + |F,(r)|2) - ^ (|G„(r)p + |F„(r)|2), 

p=l n=l 

47rrV(r) = E (|G,(r)|2 + |F,(r)|2). 

p=i 

By solving Eqs. ([8]) and ([9]) in a meshed box of size Rq self-consistently, one can calculate 
the ground state properties of a nucleus. The vector potential Vv{r) and the scalar potential 
Vs{r), energies and wave functions for bound states are also obtained. 



B. The RMF-ACCC approach 



The ACCC approach as a bound-state-type method to study resonant states in the con- 
tinuum, has the merit of numerical simplicity and easy to be implanted. The basic idea is 
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a resonant state will become a bound state if one increases the attractive potential. Res- 
onance parameters (i.e. energy and width) and wave function for a resonant state can be 
obtained by an analytic continuation from the bound-state solutions 35|. In our scheme, 
the continuation is carried out via a Pade approximant (PA). 

It should be noted that the ACCC approach is applicable for the repulsive Coulomb case, 
which has been pointed out in Ref. |35j. In this sense, it can be used for proton decay because 
of the repulsive Coulomb potential. Within the framework of the nonrelativistic Schrodinger 
equation, the ACCC approach has been applied to study single-particle resonant states for 
proton decay in spherical and deformed nuclei in Woods-Saxon potentials 



By increasing the attractive potential as Vp{r) AV^(r), a resonant state will be lowered 
and becomes a bound state if the coupling constant A is large enough. Near the branch 
point Ao, defined by the scattering threshold A;(Ao) = 0, the wave number k{X) behaves as 



k{\) 



- Ao, / > 0, 
2(A-Ao), / = 0, 



(12) 



which is an ansatz in the relativistic framework inspired by Kukiulin |35| and confirmed 



to be good by our numerical calculation 



21, 



These properties suggest an analytic 



continuation of the wave number k in the complex A plane from the bound-state region into 
the resonance region by Pade approximant of the second kind (PAII) jsSi 

k{x) ^ k^'^'^'^x) = , ^ ^ , ^ , ^T-T^' (13) 

where x = vT^^^, and co,ci. 



1 + dix + d2x'^ + . . . + d^x^ ' 
I cl, di, d2, . . . ,dN are the coefficients of PA. These co- 
efficients can be determined by a set of reference points Xi and k{xi) obtained from the 
Dirac equation with Aj > Ao, i = 1,2, ...,L + + 1. With the complex wave number 
k{X = 1) 



relation e 



resonance energy E and the width F can be extracted from the 

F 



E-t- {E,Te 



E 



and k"^ = e'^ - M^ i.e., 




/2-M, 



2 J(M2 + A;2 - kfy + Ak^k^ - 2(M2 + k^ - kj). 



In the non- relativistic limit {k <^ M), Eq. flT^ reduces to 



E 



2M 



2 

M 



(14) 



(15) 



It is evident that the continuation in the couphng constant can be replaced by the con- 
tinuation in k plane along the k{X) trajectory determined by Eq. f|T3|) to the point kn 
corresponding to the wave number for Gamow state, i.e., kji = A;[^'^1(A = 1) [35]. Similarly, 
the wave function {p{kji,r) for a resonant state can be obtained by an analytic continuation 
of the bound-state wave function ip{ki,r) in the complex k plane. One can also prove that 
the wave function ip{k,r) is an analytic function of the wave number k in the inner region 
r < Rq where the Jost function analyticity dominates [35]. Therefore, we use the technique, 
which has been adopted to find the complex resonance energy, to determine the resonance 
wave function ^pnir) = (p{kfi,r). Firstly, we construct the PA to define the resonance wave 
function at any point r in the inner region (r < Rq) 35| 

[L'N'](, , Pvik^r) ^ ao{r) + ai{r)k + a2{r)k'^ + . . . + aL'{r)k^' 

' ^ QN'ik,r)~ l + b^{r)k + b2ir)k^ + ... + bMir)k^' ' ^ ' 

where the coefficients aj(r) = 0, 1, . . . , L') and bj{r) (j = 1, 2, ... , A^') are dependent on 
r. These coefficients can be determined by a set of reference points ki and ip{ki,r) obtained 
from the Dirac equation with Aj > Aq, {i = 1, 2, L' + A^' + l). The resonance wave function 
(p{kji,r) = (p^^''^'^{k]i,r) (r < Rq) can be extrapolated in this way. It should be noted that 
the resonance wave function can also be analytically continued by the bound-state wave 
function (p{xi,r) in the complex x = y/X — Xq plane [sS^. We also try to construct the PA 
to define the resonance wave function at any point r in the inner region (r < Rq) in the 
following way, 

[L',N']( . _ PL'{x,r) _ ao{r) + ai{r)x + a2ir)x'^ + . . . + aL'ir)x^' 

The variants in above expression have the same meaning as those in Eq. (fT6|) except that 
the wave function k is replaced by x. Accordingly, the resonance wave function (pR{r) = 
(f^^' '^'^ {y/ 1 — Ao, t) (r < Rq) can be extrapolated in a similar way. Comparing the results 
of these two expressions Eq. ( |T6|) and Eq. ( |T7I) . we found they are almost the same. In the 
following calculations, the wave function will be continued by Eq. ( 1TB]1 . unless otherwise 
specified. The upper and lower components of the wave function in the inner part are both 
extracted from expression Eq. f|T6|) . 

In the outer region {r > Rq, V^^ciir) ~ ), the asymptotic wave functions of the protons 
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M 
M 





Za 
= , I.e. 

r 




Za 

r 






Za 

r 


fHt) = 





(18) 



in which a refers to fine structure constant and equals one in the natural unit {h = c = e = 1). 
The final results for the upper and lower component are respectively jsQ]. 

^ ' 2yfc^r(27+l) ^ V/ y; V/ y, / , j 

(19) 

(20) 



2v^ r(27+ 1) 



ZaSi 



K + iyM/si 



and the 



in which 7 satisfied with 7^ = — iZaY, y , 5 - 

k 'J — iy 

superscript R refers to the regular solutions of the wave function. Similarly, if 7 is substituted 
by —7, the irregular solutions of the wave function, Gj.^(A;r) and F^{kr) can be obtained 
in which IR refers to the irregular solutions. 

The outer wave function is matched to the inner wave function at r = < -Ro 

1 



y^l'^'^k^ii, r„) 



[G^iknTm) + Dikji)G'^{kRrJ] 



(21) 



where C{kji) and D{kji) = ta.n6K{kR) with ^^(^i?) the phase shift are the coefficients for 
matching. It has been found that S^ikn) is almost a constant when is large enough. Given 
the upper component, the lower component ^^(r) can be calculated from the relationship 
between the upper and the lower components derived from Eq. fllSp . i.e.. 



F.{kR,rr^) = — ^ [F^ikRT^) + D{kR)Fl''{kRrr, 
G{kR} 



(22) 



Finally, the resonance wave function is normalized according to the Zel'dovich proce- 



dures 



35| 



C. The RMF+ACCC+BCS approach 

On the basis of the RMF-ACCC approach, pairing correlation for open shell nuclei is 
treated by the BCS approximation with the width effect in resonant states of the continuum 
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considered. That is so called the RMF+ACCC+BCS approach. The BCS approximation 
has the merit of simplicity, but thought to be not reliable for those nuclei near the drip 
line, due to the improper treatment of the continuous states ^, |4o|, i.e. the inclusion of 
the spurious continuous states which change with the box size. This can be solved by the 
proper treatment of resonant continuum in the pairing correlation, taking only the real 
resonant states into consideration and kicking the spurious continuous states off from the 
continuum [?!, O]. Therefore, the contribution of the continuum, especially for those single- 
particle resonant states in the continuum, to the pairing correlation for exotic nuclei is in 
great need of consideration. 

rn 

The level density of the continuum ^41]] (single-particle resonant states) 

1 dSo, 1 r/2 



IT dSa 



TT [Br 



+ rV4 



(23) 



is introduced into the pairing gap equations instead of the discretized continuous states, in 
which 5a is the phase shift of the scattering state with angular momentum {la,ja), represent- 
ing the variation of the localization of scattering states in the energy region of a resonance. 

Suppose the pairing matrix elements are constant in the vicinity of the Fermi level 42| . 
the gap equations can be rewritten as: 



2J ^(5,-A)2 + A2 



Ei - X 



v/(5, - A)2 + A2 



2 J Jj^ v/(£„-A)2 + A- 



-dF = — 

' G 



v/(£.-A)2 + A2 



(24) 



dSry = N 



(25) 



The sums i and a run over the bound states and resonant states involved in the pairing 
calculations, and la is an energy interval related to each partial wave {la,ja)- In this way, 
the contribution of the resonant states with width effect is naturally included. For a very 
narrow resonant state, the factor becomes a delta function. 

Taking account of those resonant continuum in the gap equations, the expressions of 
various densities should be modified as: 



r)|2 - \Fa{r)\^)ga{ea)vl{ea)dea 



(26) 
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4vrrV.(r) = J](2j, + lK^(|G.(r)r + |F,(r)p) 

i 

+ J](2j, + 1) I {\G^{r)\^ + \F^{r)\^)g^{e^)vl{e^)de^ , (27) 

4vrrV3(r) = $^(2jp, + l)<(|G,,(r)p + |F,,(r)r) 
Pi 

+ 5^(2j,, + l) I {\G,Sr)? + \FUr)?)gU^r>>lS^P.)d^P. 

Pec ■^^P'^ 

-Y.('^Jn. + lK{\GnAr)\'+\FnAr)\') 

rii 

-5^(2j„^ + l) / {\G^^{r)\'+\F^^{r)\')gnAen>UenJde^a, (28) 

47rrVc(r) = $^(2j,. + l)<(|G,,(r)r + |F,,(r)r) 

Pi 

+ 5^(2j,, + 1) /" (|G,Jr)|2 + |F,,(r)n^,„(e,Jt;2je^J&,^ . (29) 

Pa •'^P^ 

The Dirac equation, meson field functions as well as the photon field function with new 
densities Eqs. f l26|27|28|l29|) are solved self-consistently in an iterative way. The total binding 
energy is given by the expression, 



E = -Enucleon + E„ + + E p + Ec + -Epair + -E'cM 



{2ji + 1) v^Ei + J2 (2ja + 1) / ga{ea)vl{ea)e^de^ 

^ J [gapsCTQ + gupcii^O + fi'pPsPo + ePc^o) d^T - j {^92^1 + ^fi'S^^O^ d^ 
^ ^ji + + + y fi'a(£:Q)t;Q(£a)Ma(£:a)(i£:a 



2 



--■41-A-1/3 
4 



(30) 



in which the last two terms are the modified pairing energy and the correction for the 
spurious center of mass motion, respectively. 

III. NUMERICAL DETAILS, RESULTS AND DISCUSSION 

Using the single-particle energies and wave functions extracted from the RMF approach 
for bound states and the RMF-ACCC approach for resonant states with relativistic Coulomb 
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wave functions boundary condition j43| , we solve the BCS gap equations with the contribu- 
tion of the continuum and obtain the Fermi energy, pairing gap as well as the occupation 
probabilities of quasi-particle states js]. The nuclear densities composed of quasi-particle 
states are recalculated with the contribution of the widths. New meson field functions as well 
as the photon field function, and new potentials are used to recalculate the single-particle 
energies and wave functions of the bound states by solving coupled Dirac equations. New 
resonance parameters and wave functions are worked out by the RMF-ACCC approach. 
Therefore, the RMF+ACCC+BCS calculations are numerically self-consistent by an itera- 
tive way until convergency. In our calculations, the energies are in the precise of 10^'^ MeV 
and that of the densities is 10~^ fm~^. 

It is our first attempt to describe the proton-rich odd nucleus in a fully self-consistent way 
by RMF+ACCC+BCS approach, with the consideration of the proton resonant continuum 
states contributions and pairing correlations. We are seriously eager to know how far we can 
go with our model. Furthermore, the recently improved experimental measurements for ^^Ne 
make it possible to give more possible comparison with our calculations. The pairing window 
is opened from all the bound states below the Fermi surface to one harmonic oscillator shell 
above the Fermi surface in our model, including the bound proton states lsi/2, IP3/25 lpi/25 
and single proton resonant states ld5/2, 2si/2. The constant pairing strength G = C/A in 
the BCS approximation, is chosen by fitting the odd-even mass difference extracted from 
three-point formula, 

Ap = ^[B{Z - 1, AT) - 2B{Z, N) + B{Z + 1, A^)], (31) 

in which the binding energies are taken from Ref. [4^]. Highly excited resonant states with 
large widths, such as ld3/2 and lf7/2 are ignored in our calculations because of the minor 
occupation probability. For odd neutrons, the blocking effect is considered in the pairing 
correlation. 

The theoretical binding energies for ^''Ne and ^^O are calculated by the RMF, RMF+BCS 
and RMF+ACCC+BCS approaches with NL3 and NLSH effective interactions 



45 



46|. The 



results of the RMF+BCS approach change with the box size because of the inclusion of 
the spurious states. Therefore, only the binding energies obtained from the RMF and the 
RMF+ACCC+BCS approaches are shown in Tab. I. The available experimental binding 
energy for ^"^Ne is 112.9276 MeV and that for ^^O is 111.955 MeV |4^. Asterisk denotes the 
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result of the RMF+ACCC+BCS calculation without resonant width contribution. In this 
case, the resonant states with widths are regarded as the discrete positive states without 
widths. It can be seen from Tab. I that the results without width effect are larger than those 
with width contribution because of the over evaluation of the pairing correlation, which is 
also mentioned in Ref. j^]- The binding energies in the RMF+ACCC+BCS calculations are 
1% larger than available data. Two-proton separation energy for ^''Ne is 0.97 MeV, and the 
corresponding theoretical one is about 0.03 MeV for NL3 and NLSH effective interactions. 
It should be noted that the binding energy for ^^O is larger than that for ^''Ne in the RMF 
calculations, which means ^^Ne is unstable against proton emission if pairing correlation is 
not taken into account. In the RMF+ACCC+BCS calculations, the binding energy for ^^O 
is less than that for ^''Ne, which means the valence protons are weakly bound in ^^Ne because 
of the attractive pairing correlation. As a consequence, pairing correlation is crucial in the 
existence of the nuclide and needs to be taken into good consideration. The contribution 
from the resonant state in the continuum to the pairing is as important as the bound states 
near the Fermi surface. The width of the resonant state has the effect of weakening the 
pairing gap and decreasing the binding energy. These factors are self-consistently included 
in our RMF+ACCC+BCS calculations so that the reasonable description of the binding 
energy can be given. 

In order to show the contribution from the single particle energies to the whole bind- 
ing energies, the single proton levels structure of ^^Ne are displayed in Tab. II, which are 
calculated by the full self-consistent RMF+ACCC+BCS approaches with NL3 and NLSH 
effective interactions, respectively. For protons, the bound states lsi/2, IP3/25 lpi/25 as well 
as the single resonant states ld5/2 and 2si/2 with narrow widths, have quite similar structure 
for NL3 and NLSH effective interactions. The single proton levels in ^''Ne are plotted in 
Fig. 1 with NL3 effective interaction to show the levels positions for each resonant states in 
potential well and barrier. The proton potential with the Coulomb barrier Vp (red curve) 
corresponds to the potential for s partial wave and Vp + l{l + l)/r^ (blue curve) refers to 
the proton potential with the Coulomb barrier and centrifugal barrier for d partial wave. It 
can be clearly seen that 2si/2 state (red line) lies just below the Coulomb barrier (/ = 0), 
and ld5/2 (blue line) lies in the Coulomb and centrifugal barrier (/ = 2). Therefore, s wave 
and d wave with positive energies lies just above the threshold and have the probabilities to 
tunnel the barrier, which are proportional to the widths of the single-particle states. 
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In Fig. 2, the density distributions in logarithm scale for the core nuclide ^^O and the two 
proton nuclide ^^Ne are presented with NL3 (solid line) and NLSH (dashed line) effective 
interactions, respectively. The error region for the experimental results are displayed by the 



gray area 



24|. 

The radii of ^^O and ^^Ne are calculated by the RMF+ACCC+BCS approach. It can 
been seen that the densities from the RMF+ACCC+BCS approach agree well with the 
experimental fit results. The matter radius of ^^Ne increases if the proton pairing correlation 
effect is taken into account. The available experimental matter radius of ^^Ne is 2.75(7) fm 
and experimental charge radius is 3.042(21) fm. From the RMF+ACCC+BCS calculations, 
the rms radius of ^'^Ne is about 2.86 fm for NL3 effective interaction, 2.78 fm for NLSH 
effective interaction; the charge radius of ^^Ne is about 3.15 fm for NL3 effective interaction, 
and 3.07 fm for NLSH effective interaction. We can conclude that the theoretical calculations 
agree well with the experimental measurements. The relative deviations are less than 5%. 

In order to investigate the detailed contribution to the radii, we plot the occupation 
probabilities of the proton single-particle states lsi/2, lp3/2, lpi/2, ld5/2 and 2si/2 in ^''Ne by 
the RMF+ACCC+BCS approach in Fig. 3. In our calculations, the occupation probability 
of (2si/2)^ is about 20%, and that for (ld5/2)^ is about 80% for the effective NL3 interaction. 



which is in accordance with the prediction in shell-model calculations |47l . |48| . The real part 
of the wave functions for the upper component of proton resonant states ld5/2 and 2si/2 
in ^^Ne by the RMF+ACCC+BCS approach are respectively plotted for NL3 and NLSH 
effective interactions in Fig. 4. The extension of the wave functions for ld5/2 and 2si/2 can 
be clearly seen at r > 3 fm, which give main contribution to the density of this range for 
i^Ne. 

Experimental studies 



22-24 



49] still can not confirm the halo structure of ^^Ne. Cor- 



responding theoretical studies are still controversial. Even if the binding energy, matter 
density and radius, agree with the experimental results, it is still hard to judge a halo 
nucleus because of small difference of the resonant levels. Therefore, the occupation prob- 
abilities of the resonant states 2si/2 and ld5/2 are discussed for the possibility of the halo 
existence. With a three-cluster model, Timofeyuk et al. obtain a (2si/2)^-dominant config- 
uration [50]. However a three-body model calculation by Garrido et al. suggested almost 
equal occupation probabilities of the (2si/2)^ and (ld5/2)^ leve^ 



mass shift, Nakamura et al. suggested a (2si/2)^ dominance 52|], but Fortune et al. oppo- 



5l| . From the Coulomb 
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sitely suggested a (ld5/2)^ dominance j47|. By calculating the interaction cross sections ai 



with a Hartree-Fock type wave function and the Glauber model [53|, Kitagawa et al. also 
proposed a (ld5/2)^ dominance. In our scheme, the occupation probability of (ld5/2)^ is 
larger than that of (2si/2)^, because the single particle level for 7r2si/2 lies just above that 
for 7rld5/2. 

For the neutron-rich nuclei close to the drip line, sd-shell structure may change so that 
the level z/2si/2 and the level z/ld5/2 inverse, which leads to the magic number 16 5J]. For 
^^Ne, if the inversion occurs, the occupation probability of 7r2si/2 will be enlarged to form 
halo. Meanwhile, the Coulomb barrier prevent the 7r2si/2 level falling. From this point 
of view, this new magic number 16 might not appear in proton-rich nuclei. At least, the 
inversion does not appear for ^''Ne in our present calculations. That means it is difficult to 
form halo. By taking some other effect into consideration, such as pNN tensor coupling 
included in the RMF calculations {ssj, the inversion of the level i^2si/2 and the level z^ld5/2 
might occur. An effort on this point within the framework of RMF theory is in the process. 



IV. SUMMARY 

The energies, widths and wave functions of the single proton resonant states for '^'^Ne are 
studied by RMF-ACCC approach with the relativistic Coulomb wave functions boundary 
condition. Pairing correlation and contribution from the single particle resonant states in the 
continuum are taken into good consideration by the resonant BCS approach, in which con- 
stant pairing strength are used. The full self-consistently microscopic RMF-I-ACCC-I-BCS 
calculations reproduce the experimental properties of two-proton halo candidate ^^Ne, such 
as the binding energies, matter radii, charge radii, densities with NL3 and NLSH effective 
interactions. The occupation probability of the proton resonant state (7r2si/2)^ is in accor- 
dance with the prediction of the shell model. The resonant states in the continuum play an 
important role in the pairing correlation as well as those bound states near Fermi surface. 
The RMF-|-ACCC-|-BCS approach can take good consideration of the pairing correlation for 
proton-rich nuclei with only a few narrow resonant states and can be safely extrapolated 
to describe exotic nuclei. Since many nuclei are deformed, it is necessary to include the 
deformation effect in the model to describe the deformed exotic nuclei. We are developing 
the deformed RMF-I-ACCC+BCS approach to study the properties of the deformed nuclei. 
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4J]. Two-proton separation 



Tab. I: Binding energies of ^^Ne calculated by the RMF and the RMF+ACCC+BCS (this work) 
approach with NL3 and NLSH effective interactions, respectively, and compared with experimental 
binding energy (S(i'^Ne)=112.928 MeV, B{^^0)=ni.955 MeV) 
energies S2p are also shown for the RMF and the RMF+ACCC+BCS approach. All energy values 
are in unit of MeV. Asterisk denotes the result of the RMF+ACCC+BCS calculations without 
resonant width contribution. 

approach -Bnls ^nl3 — -^Exp. -Bnlsh -Bnlsh - ^Exp. 

I'^Ne RMF+ACCC+BCS 114.32 1.39 113.57 0.64 

RMF+ACCC+BCS* 114.84 1.91 114.11 1.18 

RMF 112.69 -0.24 111.86 -1.07 

1^0 RMF+ACCC+BCS 114.30 2.77 113.54 2.72 

RMF 114.29 2.33 113.52 1.56 

S2p RMF+ACCC+BCS 0.02 0.03 

RMF -1.60 -1.66 
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Tab. II: Single proton energy levels structure of ^^Ne with NL3 and NLSH effective interactions, 
respectively. All energy values are in unit of MeV. 



irnlj 

E r E r 

2si/2 0.74 0.02 1.17 0.18 

ld5/2 0.52 0.002 0.66 0.0026 

lpi/2 -8.94 -9.10 

lp3/2 -15.38 -15.79 

lsi/2 -35.46 -35.98 



20 




Fig. 1: (Color online) Single proton levels in ^ Ne calculated by the RMF+ACCC+BCS approach 
for effective interaction NL3. 
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Fig. 2: Density distributions in logarithm scale for the core nuclide ^^O and ^^Ne in the 
RMF+ACCC+BCS calculations with NL3 (sohd hue) and NLSH (dashed hue) effective inter- 



actions, respectively. The gray region corresponds to the error region for experimental results 
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Fig. 3: (Color online) Occupation probabilities of the proton single-particle states lsi/2, lP3/2j 
lpi/2, Idg/a and 2si/2 in ^"^Ne by the RMF+ACCC+BCS approach. 
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Fig. 4: (Color online) Upper wave functions for proton resonant states and 2si/2 in ^^Ne by 

the RMF+ACCC+BCS approach. 
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